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Radio galaxies are among the largest and most powerful single objects known and are found at variety of redshifts, hence 
they are believed to have had a significant impact on the evolving Universe. Their relativistic jets inject considerable 
amounts of energy into the environments in which the sources reside; thus the knowledge of the fundamental properties 
(such as kinetic luminosities, lifetimes and ambient gas densities) of these sources is crucial for understanding AGN 
feedback in galaxy clusters. In this work, we explore the intrinsic and extrinsic fundamental properties of Fanaroff-Riley 
II (FR II) objects through the construction of multidimensional Monte Carlo simulations which use complete, flux limited 
radio catalogues and semi-analytical models of FR lis' time evolution to create artificial samples of radio galaxies. This 
method allows us to set better limits on the confidence intervals of the intrinsic and extrinsic fundamental parameters and 
to investigate the total energy produced and injected to the clusters' environments by populations of FR lis at various 
cosmological epochs (0.0 < z < 2.0). We find the latter estimates to be strikingly robust despite the strong degeneracy 
between the fundamental parameters - such a result points to a conclusive indicator of the scale of AGN feedback in 
clusters of galaxies. 

© .. WILEY-VCH Verkg GmbH & Co. KGaA, Weinheim 



1 Introduction 

Radio galaxies and quasars are among the largest and 
most powerful single objects known. Since they often re- 
side in groups and clusters of galaxies and are found at 
variety of redshifts these radio sources are believed to 
have had significant impact on t he evolving Univer s e (e.g . 



Gopal-Krishna & Wiita^ '200 1'; 'Lo ngair & Seldner 
Raw lings & Jarvis 2004; Zirbel 1997j). Their relativistic jets 
inject considerable amounts of energy into the ambient 
medium which surrounds active galactic nucleus (AGN); 
hence these sources are invaluable laboratories in studies of 
intergalactic and intracluste r medium, and in studies of the 
evolu tion of their hosts (e.g. lScannapieco. Silk & Bouwens 
2005h . 



Although the original models of formation and evo- 
lution of galaxy clus ters (e.g. ICowie & Binneyl 1 1977 



ixy 



Fabian & NulsenI 11977 ). which often are the habitats of 



radio galaxies, predicted self-similar growth of the struc- 
tures, it soon became apparent that they over-estimate 
cooling rates of the evolving clusters (t he cooling flow 
problem; see McNamara & NulsenI 2007 . and references 
therein); if only gravitational heating is included in the 
models, too many baryons cool to form stars and galax- 
ies become too luminous. Moreover, scaling relations 
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of galaxy clusters' properties (such as mass, tempera- 
ture. X-ray luminosity) predicted by these self-similar 
models have also been found to diverge from the ob- 
servationally derived relations (similarity breaking; e.g. 
iMarkevitchnl998 ). Mechanisms which have been proposed 
to explain these problems, referred to as non-gravitational 
heating, include supernovae and AGN feedback, sub- 
cluster merging, and thermal condu ction amo ng oth- 
ers (e.g. iMarkevitch & Vikhlinin|| 200ll; fO mma et al]i2004l; 
Ponman. Cannon & Navarroll999[Rosner & Tuckeiil989l) . 



Here we focus on the AGN feedba ck in galaxy groups and 
clusters, and on the impact of FR II dFanaroff & RilevI 19741) 



radio galaxies in particular The main effects an AGN ex- 
erts on its suiToundings include removing cold low entropy 
gas from the centres of clusters and groups to the outskirts 
of these structures. Through the work done on the environ- 
ments by the AGN jets the excess energy is transfeiTed to 
the ambient gas - this will cause an increase in the ambient 
temperature as well as in the cluster gas entropy. 

In these proceedings, we investigate the total en- 
ergy injected by FR II radio galaxies and radio-loud 
quasars into their environments . In our recent paper 



(iKapinska. Uttley & Kaiser 1 120121 hereafter KUK12) we 
presented a new numerical method which allows one to ex- 
plore the most likely fundamental parameters that define 
whole populations of FR II radio sources, i.e. distributions 
of kinetic luminosities Q, their lifetimes t,„ax and central 
gas densities po, as well as the number densities of FR II 
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radio sources, at various cosmological epochs - these have 
long been a sought-after piece of inf ormation in the study 
of radio galaxies and quasars (e.g. Barai & Wiital 12007 



DalylligQsirRawlings & Saunderslll99ll) . However. KUK12 
showed that these fundamental parameters are highly de- 
generate and to disentangle information independent con- 
straints on some of the parameters are required. These re- 
sults led to questions which we attempt to answer here: Is 
the scale of the impact of FR II radio sources on clusters ' 
environments robust despite the degeneracy? Or, can the de- 
generacies be broken ? 

The article is structured as follows. The method used 
and the results from KUK12, which this paper is based on, 
are briefly presented in fj2] The total energy produced by 
FR II radio galaxies and quasars is reported in [J3] and dis- 
cussed in [2] The summary is given in fj5] A flat Universe 
with the Hubble constant of Hq — 71 km s^^ Mpc~^, and 
r^M = 0.7 and Oa = 0.3 is assumed throughout these pro- 
ceedings. 

2 Methods 
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Fig. 1 Total energy produced by FR II radio galaxies dur- 
ing their lifetimes at different comic epochs. Data sets (con- 
sult Table [U are marked as different data points as follows: 
data set la is marked as filled circles, data set 16 as open 
circles, data set 2 as diamonds, and data set 3 as stars. For 
details see iJH 



The method which we use in our investigations is based 
on a construction of multidimensional Monte Carlo simu- 
lations, which use semi-analytical models of FR lis' time 
evolution to generate artificial samples of radio sources . 
We use the theoretical model of ^Kaiser & Alexander ( 1997 ) 
and iKaiser. Dennett- Thorpe & A lexander (1997) and wide 
ranges of searched values of the fundamental parameters 
with which the models are used in this work. Once the arti- 
ficial samples are generated one can easily find their radio 
luminosity functions (RLFs; ISchmidtill968h : these are fur- 
ther compared to the observed RLFs of FR II sources at 
various redshifts. Through the maximum likelihood method 
we found the best fitting FR lis' fundamental parame- 
ters, and their corresponding confidence intervals. We use 
complete, flux-limited radio samples, which in this work 
include t he Third Cambridge Revised Revised Catalogue 
(3CRR; Lain g, R iley & Longair 1983b arid radio sample 
constructed by |Best, Rottgerin g & LehnertI (Il999h . For all 
relevant details and assumptions adopted the reader is re- 
ferred to KUK12. 

As reported in KUK12, we found that the best fit- 
ting fundamental parameters span wide ranges of possi- 
ble values, extending over a few orders of magnitude. Al- 
though the estimated kinetic luminosities fall within the ex- 
pecte d range of 10 ^^ — lO"*^ W for FR II radio sources 



(e.g. i Rawlings & Saunders 1991), the locatio n of t he ki- 
netic luminosity break (Qb) of the Schechter ( 19761) func- 
tion (which was assumed to describe distributions of ki- 
netic luminosities) was found to be uncertain due to strong 
degeneracies between the parameters. For instance, in the 
local Universe, when allowing radio sources to have rela- 
tively short lifetimes (10^ yr) low gas density environments 
(^ 10^^^ kg m~'^) and breaks at high kinetic luminosities 
(10'^^ — 10^° W) are required in order to reconstruct the ob- 



served RLFs. Whereas for the longest lifetimes (10^ yr) of 
these sources much denser environments 10~^^ kg m^"^) 
and lower Q-q (~ 10'^* W) are necessary. Moreover, these 
parameters have also been found to evolve with cosmolog- 
ical epoch, where we assumed that the parameters change 
with redshift according to Qb(z) ~ (5b(^o)(1+^)"'' forthe 
kinetic luminosity break, and imax„, (z) = tinax^{zo){y + 
z)"' and Pin{z) = Pm(^o)(l + ^)"' for distributions of 
maximum lifetimes and central gas densities respectively, 
and where zq denotes 2 = 0. A strong degeneracy has been 
observed to occur also in nq, rit and parameters. 

Despite these degeneracies one may attempt to estimate 
the total energy (J/tot) produced by a radio source during its 
lifetime, where 



Q{t)dt. 



(1) 



In this preliminary analysis we selected four data sets to 
investigate estimates on the total energy. The data sets are 
within the 90% confidence intervals of the best fits as found 
in KUK12. We consider the best fit found (data set 16) and a 
data set in which the cosmological evolution of central gas 
densities is set to = 4.0 (discussed in KUK12, here data 
set la). In addition, we draw two data sets (2 and 3) that rep- 
resent the two extreme cases discussed above (that is cases 
when the maximum lifetime is set to the shortest and longest 
one found within 90% confidence intervals of KUK12 re- 
sults, see also Table [TJ. Furthermore, it is assumed that the 
break kinetic luminosities of each data set represent the typ- 
ical luminosities of radio sources in their respective popula- 
tions, and that radio sources at the break kinetic luminosi- 
ties have maximum lifetimes and central gas densities typ- 
ical for their populations. Based on these assumptions and 
since the kinetic luminosity is assumed to be constant over 
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Table 1 Fundamental parameters defining populations of FR II radio sources which were chosen for the study. For details 
see fj2] 



Data set 


Qb{zo) [log(W)] 
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Pm{zo) [log(kgm ^)] 




imax„(zo) [log(yr)] 
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38.0 
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0.6 


-23.3 
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-2.5 
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37.4 
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0.6 


-21.5 


5.5 


8.7 


-3.0 



the radio source's life (IKaiser & Alexander! 1 19971) . Eqn. [T] 
simplifies to 



references therein). It has been suggested before that AGNs 
with powerful radio outbursts may balance out the clus- 



Uu 



(2) 



We investigate the total power and total energy produced 
by FR lis' populations at four redshift steps, up to z = 1.5. 

3 Results 

The total energy produced by FR II radio sources for each 
data set and redshift considered is presented in Fig. [T| In- 
terestingly, the results are alike for each data set. As can 
be deduced from KUK12 results, the higher the kinetic lu- 
minosity, the shorter the lifetime of the sources is required, 
and hence, as we explicitly show here, the product of these 
degenerate results (C/tot) is robust. 

Furthermore, the total power produced by FR lis seem 
to increase with redshift. Kinetic luminosities and maxi- 
mum lifetimes of FR II radio sources undergo contrasting 
cosmological evolution, where kinetic luminosities of the 
sources increase with redshift, while their maximum life- 
times decrease with cosmic epoch. Since the cosmological 
evolution of the kinetic luminosity breaks and maximum 
lifetimes of the radio sources is known in our calculations 
(Tab. [Til the plausible redshift evolution of f/tot can be esti- 
mated. We note that Eqn. |2] can be expressed with redshift 
dependence, i.e. 



(3) 



which leads to C/tot oc (1 + z)" with n e [4.0; 6.5] in 
our current preliminary calculations. This result needs to be 
confirmed with statistics on full populations, however 

4 Discussion 

We find estimates of the total energy produced by popula- 
tions of FR II radio sources to be remarkably robust despite 
the strong degeneracy between the fundamental parameters. 
Such a result points to a compelling indicator of the scale of 
AGN feedback in groups and clusters of galaxies. 

The estimated kinetic luminosities and available energy 
of FR II radio sources suggest that radio galaxies are able to 
quench cooling fl ows at least for the duration of their life- 
times (Fig. [TJ cf. iGitti. Brighenti & McNamarall2012i and 



ters' cooling flows (e.g. Bh'zan et al.l l2004t iRaffertv et al 
2006h . The range of possible values of kinetic luminosi- 



ties (consult KUK12) may suggest that the estimates from 
our degenerate results of quenching cooling flows with ra- 
dio galaxies are less robust than those of total energy pro- 
duced (these proceedings). However, the more powerful 
outbursts are found to last for shorter time than the less pow- 
erful ones (Tab. [U. If activity in radio galaxies is recurrent, 
galaxy clusters hosting such radio sources would undergo 
alternate periods of cooling and heating, but the recurrence 
timescales may be related to the power of the radio source 
involved, where the weaker outbursts would be expected to 
occur more often. 



It is important to note here the difference between mor- 
phological classes of radio galaxies, that is FR I and FR II, 
and the possible difference in their impact on c l usters ' en- 
vironments. As pointed out by Mat hews & Guol (1201 ih Un- 
ear sizes of typical radio galaxies signific antly exceed clus - 
ters' cooling radii (rcooi ^ 100 kpc; e.g. Best et al.ll2007h . 
This is especially relevant in the case of FR II sources since 
due to their large linear extent and their strong jet head 
driven shocks one must expect most of the energy to be 
deposited beyond the cooling radius (e.g. iGitti et al.ll2007l : 
iMafliews & Guol 1201 ih . At lower redshifts (z < 1.0) it is 
the FR I sub -population that dominates investigations car- 
ried out with X-ray observations (e.g. De Young 20ll). Due 
to the difference in interaction with the ambient medium of 
the two morphological classes, it has been suggested that 
FR Is mainly contribute to heating of cluster cores, while 
FR lis contribute to overall AGN feedback, especially at 
higher redshifts where they are expected to be much more 
numerous, and may be more relevant i n the problem of 
similarity breaking (e.g. lOe Young||2010l) . Our work is cur- 
rently limited to radio sources of FR II morphology only, 
due to the availability of the theoretical models of radio 
source growth with time. In contrast to turbulent structures 
of FR Is, the idealised, and often regular, large scale ra- 
dio structures of FR lis make them relatively simple sys- 
tems to model. Although there have been attempts to de 



velop models of FR I time evolution (e.g. lBicknelllll995 



Luo & SadleJl2010l) they are still not general enough to be 



applied to whole populations of these radio sources. Once 
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such models are developed, our analysis can easily be ex- 
tended to include populations of FR I sources. 

Finally, despite the striking robustness of the total en- 
ergy estimates, one should consider possible methods for 
breaking the degeneracy between the fundamental param- 
eters; we will shortly list them here. As pointed out by 
Churazov et al.l (120021) the AGN should be finely tuned and 
neither too much nor too little heating can be accepted. 
Since the levels of required energy in order to converge ob- 
servations and theoreti cal models are determined (see e.g. 
Tozzi & Normanll200U and references therein) these esti- 



mates can be used as a tool for breaking degeneracy in the 
KUK12 results. Additionally, X-ray observations can pro- 
vide one with estimates of ambient gas densities of FR II 
radio sources at least in the local Universe, which again can 
be used to help in breaking the degeneracies. 



5 Summary 

We extend the work presented in KUK12 and investigate the 
total energy produced by populations of FR II radio sources 
during their lifetimes, and the sources' impact on their en- 
vironments, at various cosmological epochs. In these pre- 
liminary results we report that this total energy seems to 
be remarkably robust despite strong degeneracies between 
the fundamental parameters defining populations of FR II 
sources (their kinetic luminosities, lifetimes and gas den- 
sities of the environments). These results hold for all red- 
shifts considered. Moreover, we find that the total energy 
produced by FR II radio sources increase with redshift. Such 
results point to a compelling indicator of the scale of AGN 
feedback in groups and clusters of galaxies across cosmo- 
logical times. 
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